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medium.  The  mathematical  model  that  is  developed  uses  the  finite-element 
method  to  calculate  the  vibrational  characteristics  of  the  clastic  cylinder  and  the 
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SUMMARY 


PROBLEM 

Develop  a  mathematical  model  that  accurately  describes  the  electrical  and  acoustical 
response  of  a  tangentially  polarized  piezoelectric-ceramic,  staved,  free-flooded  cylinder  trans¬ 
ducer.  Implement  this  formulation  in  a  package  of  computer  programs  that,  from  a  set  of 
user-specified  dimensions  and  parameters,  will  automatically  produce  plots  of  the  complex 
electrical  impedance  and  transmitting  voltage  response  as  function;  of  frequency  and,  at  any 
selected  frequency,  will  calculate  the  vertical  directivity  pattern  and  cylinder  displacement 
distribution.  Compare  the  predicted  response  of  a  cylinder  to  that  measured  experimentally. 

RESULTS 

A  finite-element  approach  was  taken  to  the  problem  of  solving  the  equations  of 
motion  for  a  tangentially  polarized,  piezoelectric-ceramic  free-flooded  cylinder  immersed  in 
an  infinite  acoustic  fluid  medium.  The  use  of  finite-element  techniques  makes  it  possible  to 
predict  the  broadband  response  of  the  transducer.  Three  computer  programs  arc  used  in  the 
analysis:  (1 )  a  data  generation  program  which  produces  the  necessary  finite-element  input 
information;  (2)  an  instruction  program  for  the  MARTSAM  finite-clement  structural  analysis 
and  matrix  system  program;  and  (3)  a  program  which  solves  the  matrix  equations  that 
approximate  the  fluid-loaded  cylinder’s  equations  of  motion  for  the  electrical  impedance, 
transmitting  voltage  response,  and  radiated  power  as  functions  of  frequency.  In  addition, 
fora  given  frequency,  the  vertical  directivity  pattern  and  cylinder  displacement  distribution 
can  be  calculated.  Hie  ability  of  the  mathematical  model  to  accurately  predict  the  in-air 
and  in-water  responses  of  free-flooded  cylinders  was  shown  by  comparing  the  theoretical 
predictions  to  the  experimentally  measured  responses  of  three  cylinders. 
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1  INTRODUCTION 


The  free-flooded  cylinder  transducer  offers  several  distinct  advantages  over  other 
types  of  underwater  sonar  projectors.  Free-flooded  cylinders  have  good  power-to-weigbt 
ratios  and  bandwidths  and  are  also,  as  their  name  implies,  automatically  pressure  compen¬ 
sated  for  any  change  in  operating  depth.  They  arc,  therefore,  nearly  depth  insensitive  and 
have  broad  applications  whenever  a  need  arises  for  high-power  sonar  projectors  with  deep 
depth  requirements  or  for  sonar  projectors  with  very  broad  bandwidths. 

Acoustic  propagation  studies  such  as  Project  Artemis  and  LORAD  have  shown  the 
advantages  that  deep  acoustic  installations  offer  over  near-surface  locations.  Most  types  of 
transducers,  such  as  the  longitudinal  vibrator,  moving  coil,  or  piezoelectric  disk,  require 
internal  pressure  compensation  to  keep  the  exterior  ambient  pressure  from  damaging  the 
device  or  at  least  to  prevent  the  response  from  greatly  varying  with  changes  in  depth.  This 
also  provides  a  high  degree  of  desired  acoustic  decoupling  from  the  nr  lium  and  transducer 
housing  case  at  the  back  or  inside  of  the  radiating  surface.  Typical  pn  isurc  compensation 
devices  are  pressure-release  materials  such  as  air  or  other  gases,  corprti  -v  onionskin  paper, 
Min-K  and  Sonitc.  Compliant  tubes  have  also  found  great  use  as  press*.',  ivrclcase  devices, 
particularly  in  low-frequency  bender-bar  transducers.  Most  of  the**"  :*5ts,  however,  lose 
their  dynamic  compliance  under  prolonged  exposure  to  high  pre<«  .  r'  -.  depth  limited,  or 
in  the  case  of  gas  compensation,  need  complex  external  cquipirc'- 

As  one  solution  to  this  problem,  the  free-flooded  ring  <,**  <  Taler  was  developed. 
Initially,  free-flooded  cylinders  were  cast  in  one  piece  of  piczoc!  ctric  ceramic  and  were 
radially  polarized.  These  transducers  used  the  K31  coupling  of  the  electric  and  elastic  fields. 
The  use  of  vertically  striped  electrodes  made  it  possible  to  polarize  the  solid  cylinder  in  an 
approximate  alternating  tangential  polarization;  this  allowed  the  use  of  the  higher  K33 
coupling. 

When  lower  frequency  use  was  desired  and  the  larger  cylinders  became  difficult  to 
fabricate  in  a  single  piece,  a  segmented  or  staved  construction  was  designed  (Ref.  1 ).  This 
construction  allowed  the  full  use  of  the  K33  coupling.  Each  stave  was  cast  and  elcctrcdcd 
separately  and  then  glued  into  a  cylindrical  configuration.  A  drawing  of  a  typical  tangen¬ 
tially  polarized,  free-flooded  cylinder  is  shown  in  Fig.  1. 

Free-flooded  cylinders  have  been  analyzed  in  the  past  with  several  different  mathe¬ 
matical  methods  (Refs.  2-9).  These  techniques,  however,  often  could  not  predict  the  broad¬ 
band  response  of  the  cylinder  because  of  inadequacies  in  the  vibrational  models  of  the 
structure.  Finite-element  techniques  have  been  used  to  develop  very  accurate  models  for 
the  cylinder  that  can  predict  the  in-air  response  of  the  cylinder  over  a  wide  frequency  band. 
These  models  have  been  experimentally  validated  through  the  use  of  holographic  interfero¬ 
metry  (Ref.  10).  A  finite-element  model  of  a  tangentially  poiarizcd-ecramic  cylinder  has 
been  combined  with  an  acoustic  radiation  program  with  good  success  (Ref.  II).  In  this 
analysis,  finitc-clcncnt  techniques  were  used  to  calculate  the  normal  velocities  on  the 
boundaries  of  the  >  ylinder.  Then  the  Helmholtz  integral  radiation  program  CHIEF  (Ref.  1 2) 


was  used  to  calculate  tiie  acoustic  radiation  impedance  and  tltc  pressure  at  points  in  the  near 
and  far-ficlds  at  each  driving  frequency. 

Tills  last  technique,  although  successful,  is  limited  in  its  applications  due  to  the 
necessity  for  numerically  evaluating  numerous  and  expensive  integrals.  The  metliod  that  is 
devc.oped  in  this  report  uses  finite-element  techniques  as  the  mathematical  approximation 
scheme  for  the  elastic,  piezoelectric-ceramic  structure.  Finite-clement  techniques  are  also 
used  to  mathematically  model  the  acoustic  fluid  within  a  spherical  surface  which  encloses  the 
cylinder.  Boundary  conditions  which  are  simple  ratios  of  Hankel  functions  can  now  be 
applied  !c  the  spherical  surface  to  give  the  mathematical  appearance  of  a  transducer  immersed 
in  an  infinite  acoustic  fluid  medium  (frcc-flcld  response).  No  numerical  integrations  arc 
performed  for  the  acoustic  radiation  part  of  the  problem  since  the  wave  equation  separates 
on  a  spherical  surface  when  the  pressures  and  normal  velocities  am  expanded  in  spherical 
harmonic  scries. 

Section  2  of  this  report  presents  the  derivation  of  the  matrix  equations  that 
represent  the  equations  of  motion  of  a  free-flooded  cylinder.  Techniques  for  calculating 
both  the  cylinder’s  in-air  and  in-water  responses  arc  given. 

Section  3  describes  the  computer  programs  that  have  been  developed  to  implement 
the  mathematical  techniques  described  in  Section  2.  Front  a  set  of  user-specified  dimensions 
and  parameters,  three  programs  automatically  generate  the  necessary  flnitc-dcmcnt  idealiza¬ 
tions  for  the  cylinder  and  surrounding  sphere  of  acoustic  fluid,  set  up  (lie  required  matrix 
equations,  and  then  calculate  the  complex  electrical  impedance,  transmitting  voltage  response 
and  radiated  power  as  functions  of  frequency.  In  addition,  at  any  frequency',  the  vertical 
directivity  pattern  and  the  displacement  distribution  on  the  cylinder  can  be  calculated. 

Section  4  compares  the  theoretically  calculated  and  experimentally  measured 
responses  of  three  free-flooded  cylinder  transducers.  Comparisons  are  made  between  the 
in-air  resonance  and  antircsnnancc  frequencies,  electrical  capacitance,  in-water  electrical 
impedance,  transmitting  voltage  response,  and  directivity  pattern*.  Predictions  arc  nude 
for  maximum  source  level  and  maximum  radiated  power  at  selected  frequencies.  The  diffi¬ 
culties  involved  in  predicting  the  response  of  a  ceramic  cylinder  housed  in  a  caster-ofi-filled 
rubber  and  brass  vase  by  means  of  a  mathematical  model  for  only  the  bare  ceramic  cylinder 
arc  discussed. 


2  FORMULATION 

All  of  Euclidean  3-spa cc  is  divided  into  two  regions.  R|  and  R->  (Fig.  2).  Region  R| 
has  a  finite  voiumc  and  contains  all  sources  of  sound;  i.c..  R  j  contains  the  free-flooded 
cylinder  and  a  portion  of  the  acoustic  fluid.  This  region  is  bounded  by  the  spherical  surface 
S] .  Region  R?  is  infinite  in  extent  and  b  filled  with  a  homogeneous  acoustic  fluid.  The 
outward-drawn  normal  from  a  point  {  on  surface  Sj  b  defined  as  n*  The  vector  d(x4)  is 
the  distance  from  a  point  {  on  Sj  to  the  field  point  x  in  region  Rj.  Tire  point  x  b  interior 
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to  S|  in  region  Rj .  Because  the  surface  S j  is  spherical,  simple  analytical  expressions  can  be 
found  for  the  relationship  between  the  surface  pressures  and  normal  velocities.  This 
acoustic  radiation  impedance  is  the  bouisdary  condition  that  must  be  applied  to  the  equa¬ 
tions  of  region  R  j  to  yield  the  interaction  of  region  R  j  with  region  Rj.  Hence,  the  solution 
of  the  coupled  struct  ure-vibrat  ion/acouslic-radiation  problem  can  be  reduced  to  a  system  of 
boundary-vjlue  equation*  whose  spatial  domain  is  now  a  finite  volume.  A  more  detailed 
explanation  of  this  formulation  can  be  found  in  Ref.  13.  An  approximation  method  is 
now  used  to  solve  the  equations  of  region  Rj  subject  to  the  analytical  boundary  conditions 
for  surface  S|.  This  approximation  scheme  utilizes  the  finite-element  method  (Refs.  14-16). 
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In  the  finite-element  formulation,  tile  equations  of  motion  for  an  elastic  solid  can 
be  written  in  the  matrix  form 


|MSJ  {U}+|KSI  {u}=  {U }  , 


(1) 


where  (MSI  is  the  structural  mass  matrix,  IKj)  the  structural  stiffness  matrix,  {U}  the 
generalized  nodal-point  forces,  and  {u}  the  nodal-point  displacements.  The  derivative  0  is 


32u 


defined  as  — r- .  Tlie  structural  mass  matrix  is  defined  by 
3t2 


IMJ  2  /  /{N'1T(N'I  dV  „ 


(2) 


where  p  is  the  density  of  the  piezoelectric  ceramic  and  IN')  the  structural-element  inter¬ 
polation  function.  Hie  structural  stiffness  matrix  is  defined  by 


IKJ  =  /ib')TID)  IB')  dV  , 


(3) 


where  IB' I  is  the  matrix  of  derivatives  of  the  structural-element  interpolation  function 
(the  matrix  that  relates  the  strains  to  the  nodal-point  displacements),  and  ID)  denotes  the 
elasticity  matrix  containing  the  appropriate  material  properties. 


The  generalized  nodal-point  forces  can  be  written  as  the  sum  of  two  forces. 


(U}=  {P}  +  {FS}  . 


(4) 


where  (P)  is  a  vector  consisting  of  the  nodal-point  forces  arising  from  the  fluid  pressure  on 
the  fluid-structure  interfaces,  and  (F$  }is  a  vector  consisting  of  externally  applied  nodal- 
poir.t  forces  (in  this  ease,  an  electrical  forcing  function  applied  to  the  piezoelectric  structure 
via  a  sinusoidally  varying  charge). 


Similarly,  the  equations  of  motion  which  govern  the  behavior  of  the  acoustic  fluid 
in  region  Rj  are  formulated  in  Refs.  17  and  18  as 


IMf|  {p}  +  lKfI  {p}={C}  , 


(5) 


where  |Mf)  is  the  fluid  mass  matrix,  (Kf)  the  fluid  stiffness  matrix.  (C)  the  generalized 
nodal-point  forces,  and  (p)  the  nodal-point  pressures.  Again,  p  is  defined  as  32p  .  The 
fluid  mass  matrix  is  defined  by  31 2 


?MfJ  4/lNITIN)  dV  , 


(6) 


a 


gl 
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where  c  is  the  speed  of  sound  in  the  fluid  and  [N]  the  fluid-element  interpolation  function. 
The  fluid  stiffness  matrix  is  defined  by 


lKf] 


/[B1T[BJ  dV  , 
V 


(7) 


where  [B]  is  the  matrix  of  derivatives  of  the  fluid-element  interpolation  function. 

The  generalized  acoustic  fluid  nodal-point  forces  can  also  be  written  as  the  sum  of 
two  forces. 


{C}  =  {S}+  {Ff }  ,  (8) 

where  {S}  is  a  vector  consisting  of  the  nodal-point  forces  due  to  the  acceleration  of  the 
boundaries  of  the  clastic  structure,  and  {Ff }  is  a  vector  consisting  of  the  nodal-point  forces 
that  arise  from  the  interaction  between  the  acoustic  fluid  in  region  Rj  and  the  acoustic 
fluid  in  the  infinite  exterior  region  R2. 

If  the  pressure  interpolation  function  Nj(x,y,z)  for  the  acoustic  fluid  elements  is 
defined  as 


p(x.y,z)  =  2  Nj  (x.y,z)  pj 


(9) 


where  pj  is  the  pressure  at  the  i1*1  node  and  p(x,yTz)  is  the  pressure  anywhere  in  rhe 
clement,  then  Ref.  17  defines  the  fluid  forces  at  node  i  as 


elements 


(10) 


3P  . 


where  — —is  the  normal  derivative  of  the 
«n 


pressure.  Using  the  boundary  condition 


9P 

dn 


=  -pv" 


where  vn  is  the  norma!  velocity  of  the  fluid,  Eq.  (10)  can  be  written  as 


*s>i  =  -P  Z 

elements 


« 


(11) 


(12) 


8 


But  vn  itself  can  be  written  in  terms  of  the  nodal-point  accelerations  as 


vn  =  ?  Nj  ii:  , 


‘J  J 


(13) 


where  Nj  is  the  anpropriate  structural-element  interpolation  function  for  the  normal  dis¬ 
placement  at  the  structure-fluid  boundary,  and  ii:  is  the  acceleration  of  the  j^1  nodal  point 
of  the  structure.  Combining  Eqs.  (12)  and  (13),  (S;)  now  has  u  e  form 

ip  1  — % 


Wi'-ol  s  /NiNjfljdS  . 

dements  j  Lb  J 


(14) 


Similarly,  the  fluid -stiueture  forcing  function  in  Eq.  (4)  can  be  written  as 


:T">. 


Z  s  f  Nj  Nj  Pj  dS 
elements  j  L  S 


(15) 


The  force  ierms  (S)j  and  (F)j  are  the  coupling  between  the  acoustic  fluid  and  the 
elasti;  straw  ure.  Comparing  Eqs  (14)  and  (15)  reveals  that  by  defining  a  coupling  matrix 
[V]  such  that 


<T),j  Z  /Ni  Nj  dSl  • 
,J  elements  LS  J  J 


(16) 


Eqs.  (14)  and  (1 5)  can  be  written  in  matrix  form  as 
{S}=  -/>fT]T{u}  , 
and 

(P)  =  IT]  {p}  . 


(17) 


(18) 


The  combined  structure-vibration/acoustic-radiation  problem  can  now  be  written 
as  two  sets  of  coupled  equations.  The  structure  equations  of  motion  appear  as 


IMS]  {«}  +  [Ks]  {u}  =  [T]  CP)  +  (Fs)  , 
and  the  fluid  equations  appear  in  the  form 

[tf/1  W  +  [Kf]  {p}  =  -p[T]T  (ii)  +  {Ff}  . 


(19) 


(20) 


IN-AIR  STRUCTURAL  RESPONSE 


The  structure  being  considered  here  is  composed  of  a  material  having  not  only  elastic 
properties  but  piezoelectric  properties  as  well.  Reference  1 9  formulates  the  equations  of 
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motion  for  a  piezoelectric  solid.  The  assembled  set  of  structural  equations  can  be  written 
in  the  form 


Muu  0 


where  {v}  is  the  vector  of  nodal-point  electrical  potentials,  {Q }  the  vector  of  electrical 
charges,  [Kuu]  replaces  [Ks],  [Muu]  replaces  [Ms],  [K^J  is  the  piezoelectric  stiffness 
matrix,  and  [K^]  and  [K^u]  are  the  electro-elastic  coupling  matrices. 

Equation  (21)  can  now  te  used  to  predict  the  in-air  response  of  the  free-flooded  cylinder 
transducer  by  calculating  the  resonance  (short  circuit)  and  antiresonance  (open  circuit) 
frequencies  and  electrical  capacitance. 

Ihe  electrical-potential  degrees  of  freedom  must  first  be  condensed  to  one  degree  of 
freedom  by  applying  the  appropriate  boundary  conditions.  This  one  electrical-potential 
degree  of  freedom  is  the  potential  which  is  applied  across  the  transducer’s  input  terminals. 

The  potential  degrees  of  freedom  for  the  nodal  points  on  one  side  of  a  stave  of  the  cylinder 
are  all  set  equal  to  the  driving  potential  Vq,  and  the  potential  degrees  of  freedom  on  the 
opposite  side  of  the  stave  are  all  set  to  zero.  These  two  sides  are  the  silvered  surfaces  (see 
Fig.  1)  to  which  the  driven  and  grounded  input  terminals  are  attached.  The  remaining 
potential  degrees  of  freedom  in  the  stave’s  interior  are  left  free.  The  application  of  these 
boundary  conditions  and  a  repartitioning  makes  it  possible  to  write  Eq.  (21 )  in  the  form 
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The  second  equation  of  Eq.  (22)  can  be  written  as 

1*2,1 1  "o  +  1*2, 2l  W1  • 


w  =  -  (K2i2r>  ik2i,i  “o  . 

Substituting  Eq.  (23)  into  the  first  equation  of  Eq.  (22)  yields 


id  MK,,,]  L,  -lK!>2HK2j2rMK2>1] 


0 

Q  > 
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Again,  assuming  a  harmonic  time  dependence,  Eq.  (29)  becomes 
f-“A  Muu  +  Kuu]  {uA}  =  {0}  , 


(31) 


and  the  antiresonance  frequencies  and  their  associated  modeshapes  can  now  be  found.  This 
is  the  open-circuit  case. 


For  a  free-flooded  cylinder,  three  resonant  modes  of  vibration  are  of  interest:  the 
breathing,  the  bending,  and  the  length  modes.  Fig.  3  shows  these  three  modes  in  a  cross- 
sectional  display  of  the  cylinder.  The  breathing  mode  is  the  best  radiator  of  sound  and  is 
strongly  excited  in  air.  The  bending  mode  is  weakly  coupled  to  the  electrical  drive  and  is 
weakly  excited  in  air.  The  length  mode  is  a  poor  radiator  of  sound  because  of  the  small 
volume  displacements  involved.  When  the  cylinder  is  placed  in  water,  however,  combinations 
of  the  free-vibrational  modes  of  the  cylinder  may  become  excited  because  of  the  acoustic 
ladiation  coupling. 


BREATHING 


Figure  3.  Cross-scclional  drawing  of  the  first  three  resonant  modes  of  a  frc-r-flooded  cylinder. 
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The  electrical  capacitance  at  a  circular  frequency  a>c  can  be  calculated  by  rewriting 
Eq.  (26)  in  the  form 


^uu  *  wc  ^uu 


\u  =  0 

Kw  I  **0  Q 

J  %  '  ' 


Solving  the  first  equation  of  Eq.  (32)  for  {u}  yields 


{u}=-tKuu-<o2Muur1  (Kyy)  vq  . 

Substituting  Eq.  .33)  into  the  latest  equation  of  Eq.  (32)  results  in 

(  (v-n,>  IKuu-“c  Muul_l  (*„,>♦  Kw)  -0  =  Q  , 


j 


Cvo  =  Q  > 

where  the  capacitance  C  is  defined  by 

C  =  K„„  -  <Km)  Muur‘  {Km)  .  (35) 

IN-WATER  RESPONSE 

If  the  in-air  results  calculated  by  means  of  Eqs.  (27),  (31 )  and  (35)  are  as  desired, 
then  this  finite-element  mathematical  model  for  the  free-flooded  cylinder  transducer  can 
now  be  combined  in  a  mathematical  model  of  the  transducer  immersed  in  an  infinite  fluid 
medium;  i.e.,  the  in-water,  frce-field  response  of  the  transducer  can  now  be  calculated. 

Equations  (  20)  and  (26)  can  be  combined  into  one  matrix  equation  of  the  form 

*Mf  pTT  0*1  p  1  TKf  0  0  I  p  Ff 

0  Muu  0  ti  +  -T  Kuu  Ku„  u  =  0  .  (36) 

.0  0  oj[%J  Lo  Km  kJU  |q 

The  nodal-point  forces  Fj- acting  on  the  surface  Sj  arise  from  the  presence  of  the 
exterior  infinite  acoustic  fluid  of  region  R2.  To  obtain  an  expression  for  this  forcing 
function,  the  surface  Helmholtz  integral  equation  can  be  used  (Refs.  20-22).  This  can  be 
considered  as  the  boundary  condition  that  is  applied  to  surface  S  j  because  of  the  interaction 
between  region  R j  and  region  R2.  The  integral  formulation  relates  the  acoustic  pressure  p 
at  any  point  x  in  the  region  R2  to  the  integral  of  the  acoustic  pressure  and  its  normal 
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derivative  o\er  the  closed  surface  Sj.  Th'  surface  Helmholtz  integral  equation  can  be 
written  in  the  form 


1  f  (  3  re-jkd(x;£)l 

4jt  s,  l  3^  [  d(x,£)  J 

re-jkd(x,i)‘l  . 

poarj } dS  (i)  • 


where  k  is  the  wave  number  cj/c  (co  the  circular  frequency  and  c  the  speed  of  sound  in  the 
fluid),  and  the  pressure  in  region  R->  must  satisfy  the  radiation  condition  on  a  spherical 
surface  Sr  a?  infinity  (Ref.  21 );  i.c., 

/  m  -  jkp<r)  2  „s  =  0  .  (38) 

R->o>  J  3r 

SR  r=R 

Tlu  relationship  between  the  pressure  field  and  its  normal  derivative  on  the  surface 
S|  can  be  calculated  in  several  ways  (Refs.  1 2,  23  and  24).  However,  since  the  surface  S| 
has  been  chosen  to  be  spherical,  the  steady-state  acoustic-radiation  Green’s  function 
separates  v'iicn  expanded  in  spherical  harmonics  (Ref.  25),  and  the  radiation  impedance 
associated  vith  each  spherical  harmonic  subspacc  can  be  obtained  in  an  analytical  form. 
This  avoid?'  the  costly  numerical  evaluation  of  numerous  integrals,  such  as  is  performed  in 
the  formulation  given  in  Ref.  1 1 .  By  expanding  the  pressures  and  normal  velocities  at  the 
fmite-clen.ent  nodal  points  on  surface  Sj  in  a  finite  spherical  harmonic  series,  it  is  possible 
to  derive  an  analytical  solution  for  the  radiation  impedance  on  S  j .  Reference  1 3  gives  the 
details  of  this  calculation.  Writing  Eq.  (1 1 )  as 

=  -j«pvn(£)  ,  (39) 


the  relationship  between  the  coefficients  of  a  spherical  harmonic  expansion  of  the  pressure 
and  the  coefficients  of  normal  velocity  on  the  spherical  surface  S  j  of  radius  a  can  be 
written  as 

v£  (a)  =  t  *  r  —  ;  Pp  (a)  -  (40) 

*  jtopa  A^tka)  K 

The  radiation  admittance  1  A^(ka)  is  defined  in  Appendix  A  of  Ref.  13  as 


IifcVka) 

l,X*(ka>5  i+,-ka,^ 


where  h^(ka)  is  the  £Vl  order  spherical  Hankcl  function  (Ref.  26).  The  expansion 
coefficients  p^(a)  and  v”(a)  arc  the  £^1  components  of  the  surface  pressure  and  normal 
velocity,  respectively,  and  are  given  by 

oo 

P(a,O,0)  =  £  p£(a)Y0(0^)  ,  (42) 

£= 0 


and 


with 


vn(a ,0,0)  =£  v"  (a)Y<>  (0,0)  , 


1=0 

2i r  jt 

/ 

0  0 


P^a)  =  /“*/  p(a,0,0)Y®  (0,^)  sin(0)d0  , 


(43) 


(44) 


and 


2ir  ,  N 

."(a)  =fd<pf  vn  (a ,0,0)  Y°  (0J)  sin(0)  d6  , 


(45) 


where  (a, 0,0)  arc  the  spherical  coordinates  of  the  point  £  on  S  j . 

In  practice,  the  summations  in  Eqs.  (42)  and  (43)  arc  carried  out  over  only  the  first 
few  terms.  Pic  exact  number  of  terms  th3t  is  necessary  can  be  found  by  first  calculating 
the  maximum  value  of  ka  and  then  finding  the  maximum  value  of  £  for  which  the  imaginary 
part  of  the  acoustic  admittance  is  zero  below  this  ka  value.  Figure  4  is  a  graph  of  the  greatest 
value  of  ka  versus  the  greatest  value  of  £.  For  example,  for  a  cylinder  to  be  operated  to  a 
maximum  frequency  of  20  kHz  and  enclosed  in  a  fluid  sphere  0.1  m  in  radius,  the  greatest 
2rrfa 

value  of  ka  =  — ™  is  less  than  ten.  From  Fig.  4,  the  greatest  value  of  £  (denoted  bv-C  )  is 

therefore  fourteen.  However,  only  the  even-order  terms  need  be  retained  because  the 
orientation  and  electric  drive  configuration  of  the  cylinder  produce  an  even  symmetry  plane 
at  0  =  90  deg.  Tliercforc,  only  those  spherical  harmonics  which  are  even  functions  under  a 
reflection  through  this  symmetry  plane  arc  retained.  Tire  summations  for  the  above  example 
contain  only  eight  terms. 

Pic  form  of  the  forces  acting  on  the  fluid-sphere  surface  Sj  is  the  same  as  that  of 
the  for  *cs  due  to  the  presence  of  the  vibrating  cylinder  that  act  on  the  acoustic  fluM  }Eq.  (1 7)  J 
except  that  the  coupling  matrix  now  couples  the  infinite  acoustic  fluid  exterior  region  to 
the  acoustic  fluid  sphere  via  the  surface-fluid  finite-element  interfaces.  Ptis  coupling  matrix 
will  be  denoted  by  (Ms)  since  its  elements  arc  analogous  to  a  consistent  mass  matrix;  c.g., 
they  arc  the  integrals  of  a  quadratic  form  of  the  interpolation  function  Nj  over  the  closed 


where 


{Fs}  =  -plMjjj  {v11}  .  ,  (48) 

•  3vn 

The  vector  {v11}  contains  the  accelerations  defined  by  for  each  nodal  point  on  Sj. 

Tlie  subset  of  spherical  harmonics  that  has  been  chosen  must  now  be  normalized  to 
be  compatible  with  the  finite-element  idealization  of  surface  S| .  The  normalized  spherical 
harmonics  that  arc  the  elements  of  the  linear  transformation  n.  .trix  [  L)  arc  given  by 


W  =  ^  v?  <*i>  • 


(49) 


and  the  full  transformation  matrix  appears  in  the  form 

Lq  (0)  L2{0)  L4(C)  -  -  -  Lx  (0) 

I  LI  =  I  (Oj)  lo  (0t)  L4 15;)  -  •  L t  (Oj) 

•  -  •  * 

Lfl  (x/2)  1-2  (*/2)  L4  (»/2) '  -  L£  0/2) 

Hie  transformation  matrix  has  been  normalized  such  that 

|L)T|MS)[L!  =  U1  , 

where  [I]  is  the  identity  matrix. 


(50) 


*51) 


To  transform  the  equations  for  surface  Sj  to  the  chosen  subset  of  normalized 
spherical  harmonic  basis  vectors,  the  right-land  side  of  Eq.  (47)  is  multiplied  by  a  transfor¬ 
mation  matrix  such  that 


(52) 


The  spherical  harmonic  expansion  of  the  derivatives  of  the  nodal-point  velocities  can  be 
written  in  the  matrix  form 
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Using  Eqs.  (5 1 )  and  (53).  Eq.  (52)  becomes 
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(54) 


Substituting  Eq.  (40)  for  the  normal  velocity  components  in  Eq.  (54)  yields 
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(55) 


Similarly,  the  left-hand  side  of  Eq.  (47>  can  be  multiplied  by  the  transformation 
matrix  and  the  surface  pressures  expanded  in  spherical  harmonics. 


such  that 


(p}s-  I  LI  (P£> 
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Combining  Eqs.  (55)  and  (57),  the  transformation  of  Eq.  (47)  can  now  be  expressed  as 
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.  (58) 


Equation  (58)  can  be  solved  at  a  set  of  circular  driving  frequencies  ajfor  the  unknown 
Held  variables  p? .  v.  p-  and  u  in  terms  of  the  known  electrical  forcing  function  Q.  However, 
this  is  a  large  system  of  equations,  and  a  total  solution  at  each  frequency  would  be  costly. 
Typically,  the  equations  arc  more  than  one  hundred  fifty  in  number.  To  reduce  the  number 
of  equations  that  need  to  be  solved  at  each  frequency,  Eq.  (58)  is  partitioned  as  shown  above. 
Tiic  second  set  of  equations  of  Eq.  (58)  can  be  solved  for 


and  the  solution  substituted  into  the  first  set  of  equations.  Partitioning  the  matrices  such 
tliat  the  interior  problem  (Refs.  1 2  and  20)  can  be  solved  separately  yields  a  small  system 
of  surface  equations  to  be  solved  at  each  driving  frequency.  If  the  driving  frequency  is 
near  an  cigcnfrequcncy  of  the  interior  problem,  then  the  surface  equations  are  overdetermined 
by  one  degree  of  freedom.  That  one  degree  of  freedom  is  the  eigenvector  of  the  interior 
problem  that  is  associated  with  tliat  eigen  frequency. 


To  calculate  the  interior  eigenvalues  and  eigenvectors,  the  lower  right-hand  blocks  of 
Eq.  (58)  can  be  written  separately.  However,  a  simple  modification,  suggested  by  O.  C. 
Zicnkiewicz  and  B.  M.  Irons  (Ref.  1 7).  produces  a  symmetric  form  for  the  equations 
describing  the  interior  problem  which  makes  it  caster  and  less  expensive  to  calculate  all  of 
the  eigenvalues  and  eigenvectors.  Solving  the  third  equation  of  Eq.  (58)  for  {p  }  j  yields 


{pjj  -  -IKrl-l  IMflu  (LI  (Pi)  -  (Kfljj  IMflj,  {p)j 
-IKfljj  IpT|T  <U>  -  [Kflj-j  [KfJjj  ILI  (Pj) 


(59) 
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(67) 
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Equation  (66)  can  be  solved  in  the  frequency  domain  by  first  solving  the  second 
set  of  equations  for  {g}  and  substituting  the  solution  into  the  first  set  of  equations.  For 
driving  frequencies  that  are  close  to  an  interior  eigenfreauency,  a  solution  to  the  first  set  of 
equations  cannot  be  found  because  factors  of  l/(u>?  -  are  present.  To  avoid  this  singularity 
problem,  the  transformed  surface  ptessure  and  electric  potential  equations  arc  augmented  by 
the  eigenvector  that  is  associated  with  the  interior  cigenfrcquency  that  is  close  to  the  driving 
frequency.  Therefore,  the  system  of  equations  never  needs  to  be  overdetermined  by  more 
than  one  equation.  However,  for  convenience,  all  those  equations  that  are  associated  with 
eigenfrequencies  of  the  interior  problem  that  lie  inside  the  frequency  band  of  interest  can 
as  a  block  be  used  to  augment  the  surface  equations;  i.c.,  [Kj  j  1  and  [Mj  j  J  are  enlarged 
by  k  rows  and  columns.  Equation  (66)  now  appears  as 
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if  the  following  definitions  arc  made 
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,  ,  =  (Ki.2>n-k 

1,2  (_  a 


Ik2,)]  =  «K2,lVk  01 


The  subscripts  k  and  n-k  refer  to  the  blocks  of  the  matrices  that  are  associated  with  the  k 
smallest  and  n-k  largest  (the  remainder  out  of  n  total)  interior  eigenvalues.  The  value  of  k 
is  associated  with  the  lowest  eigen  frequency  that  lies  above  the  frequency  band  of  interest. 


Assuming  a  time  dependence  of  eJU)t,  the  second  equation  of  Eq.  (68)  can  be  solved  for 


?tl‘k}  [wn-k  -  w2  ] 


[o)2M2>i  -  K2>1  ]  v 
.  8k 


Substituting  Eq.  (70)  into  the  solution  to  the  first  equation  of  Eq.  (68)  results  in 


_  _  _  /  1  \ 

v  =  f(-to2 Mj (1  +  K1>1)+  (-w2  Mj  2  +  Kj  2)  f  j  2  ) 
gi,  X^n-k  w  /Dia 


•  (u)2  M21  -  K2  I)J'l  |  Q 


In  general,  this  set  of  equations  is  not  large;  therefore,  solving  the  system  at  a  great  number 
of  frequencies  is  not  expensive. 


The  complex  electrical  impedance  of  the  piezoelectric-ceramic  cylinder  can  now  be 
calculated  as  a  function  of  frequency.  Defining  the  complex  impedance  Z  from 


i>  =  ZI  , 


where  the  current  I  is  defined  by 


3t  ’ 


3  1 


and 


Q  =  Q0eja)t  , 

the  electrical-charge  forcing  function  can  be  written  as 


Q  =  -A-  • 

JO) 


(73) 


Substituting  Eqs.  (72)  and  (73)  into  Eq.  (71 )  yields 


« 

% 


m 


P£ 

Z 

gk 


-  ((-a,2  M t >t  +  K,  _ t )  +  (-w2  M !  >2  +  K i  >2) 


Diag 


.(u2m2)I-k2>1)]-1 


0 

1/jw 

0 


(74) 


The  solution  vector  for  Eq.  (71 )  that  is  obtained  for  constant-voltage  drive 


P£ 
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can  be  calculated  from  Eq.  (74)  by  normalizing  the  solution  to  Eq.  (74)  to  the  impedance 
Z,  i.e., 
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(75) 


Hie  surface  and  interior  acoustic-fluid-sphere  pressures  and  the  cylinder  displacements 
per  input  volt  can  now  be  calculated  from  Eqs.  (56),  (65)  and  (70). 

Since  the  mathematical  formulation  presented  here  includes  no  structural  damping, 
the  efficiency  of  the  txarsducer  is  100  percent.  Therefore,  the  radiated  acoustic  power  is 
equal  to  the  electric  power  delivered  to  the  transducer.  The  radiated  power  P  can  be  calculated 
from 


P  =  Vi  Re  (*>!*). 


(76) 
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Substituting  Eq.  (72)  into  E}.  (76)  yields 

P  =  Vz  Re  ( vv*Z *)  .  (77) 

When  vv*  =  1 ,  Eq.  (77)  calculates  the  radiated  power  per  input  volt  and  can  be  written  as 
Re(Z) 


P  =  Vz 


[Zj- 


The  components  of  a  spherical  harmonic  expansion  of  the  pressure  at  any  point 
x(r,0')  in  region  R2  are  given  by 

,  ,  ,  .  h^(kr) 

^(r)  =  P^h|^  ’ 


(79) 


where  p^(a)  is  given  by  Eq.  (75).  The  actual  pressure  at  x  can  now  be  calculated  by  using 
Eq.  (56)  for  the  appropriate  angular  coordinates;  i.e., 


{p(r,0')l  =  [L(0'11  (p^(r)}  . 


(80) 


The  derivations  of  these  equations  are  found  in  Appendix  B  of  Ref.  13.  If  the  transforma¬ 
tion  matrix  f  L]  of  Eq.  (56)  is  now  used  in  Eq.  (80),  then  the  pressures  in  region  R2  are  cal¬ 
culated  at  points  on  the  surface  of  a  sphere  of  radius  r  with  the  same  angular  distribution  as 
the  nodal  points  on  surface  S  j .  The  calculated  vector  {p(r,0')}  is  really  a  vertical  directivity 
pattern. 


3  COMPUTER  IMPLEMENTATION 

The  mathematical  formulation  that  has  been  constructed  for  the  analysis  of 
tangentially  polarized,  piezoelectric-ceramic,  free-flooded  cylinder  transducers  was 
implemented  ;:s  a  series  of  three  computer  programs.  The  first  program  generates  input  data 
for  the  second  and  third  programs.  The  second  program  is  a  MARTSAM  finite-element 
instruction  program,  and  it  calculates  the  cylinder’s  in-air  response  and  forms  the  matrices 
that  comprise  Eq.  (74).  The  third  program  performs  a  frequency  sweep,  and  at  each  frequency 
it  calculates  the  complex  electrical  impedance,  radiated  power  per  input  volt,  and  transmit¬ 
ting  voltage  response.  Directivity  patterns  and  cylinder  displacement  distributions  can  be 
calculated  for  selected  frequencies. 
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PROGRAM  DATFFC* 


This  program  takes  given  dimensions  of  the  cylinder  to  be  analyzed  and  specified 
finite-element  idealizations  of  the  acoustic-fluid-sphere  surface  and  the  cylinder  and  generates 
the  necessary  data  for  the  second  and  third  programs.  The  computer  program  DATFFC 
needs  to  have  the  following  input  parameters  specified  by  the  user  in  MKS  units: 

RI,  inner  radius  of  the  cylinder 

RO,  outer  radius  of  the  cylinder 

RLNTH,  total  length  of  the  cylinder 

NSTAVE,  number  of  staves  that  comprise  the  cylinder 

SMALLA,  radius  of  the  fluid-sphere  surface  Sj 

RHO,  density  of  the  acoustic  fluid 

RCV,  speed  of  sound  in  the  acoustic  fluid 

DENSTY,  density  of  the  piezoelectric-ceramic  material 

RMATPM(I),  ten  piezoelectric-ceramic  material  parameters  (see  Ref.  8) 

NPSR,  number  of  fluid-sphere  surface  elementr 
NLR,  number  of  structure  elements  along  the  half-length  of  the  cylinder 
NTR,  number  of  structure  elements  through  the  thickness  of  the  cylinder 
NSPHH,  number  of  terms  used  in  the  spherical  harmonic  expansions 
KAUG,  number  of  degrees  of  freedom  used  to  augment  the  surface 

equations  to  prevent  singularities  in  the  solution  over  the  frequency 
band  of  interest 

CAPFRQ,  the  frequency  at  which  the  in-air  electrical  capacitance  of  the 
transducer  is  to  be  calculated 

The  structure  has  midplane  reflective  symmetry  in  the  z-direction  and  N-fold 
reflective  and  rotational  symmetry  about  the  z-axis.  Tire  structure  therefore  consists  of  N 
unit  cells.  It  is  necessary  to  mathematically  model  only  one  unit  cell  and  apply  the  boundary 
conditions  of  zero  azimuthal  displacement  for  the  boundaries  of  the  unit  cell  that  are 
shared  by  its  neighboring  cells.  Because  the  electrical  driving  force  is  axisymmetric,  only 
m  =  0  circular  harmonic  motion  is  considered  (see  Ref.  1 3).  This  means  that  the  cylinder 
displacements  and  the  acoustic-fluid  pressures  arc  independent  of  the  azimuthal  angle  <(>,  and 
if  the  displacements  of  the  cylinder  are  written  in  terms  of  a  cylindrical  coordinate  system, 
then  all  azimuthal  degrees  of  freedom  can  be  set  to  zero  with  no  loss  of  accuracy  over  the 
frequency  band  of  interest.  The  mathematical  modeling  can  therefore  be  reduced  to  a  two- 
dimensional  problem  even  though  the  unit  cell  of  the  cylinder,  one  stave,  is  not  axisymmetric. 
In  addition,  because  the  cylinder  has  midplane  reflective  symmetry  (see  Fig.  5),  modeling  of 
only  the  top  half  of  a  stave  and  the  top  half  of  the  fluid  sphere  and  applying  the  symmetry 
boundary  condition  of  setting  to  zero  the  z-displacements  of  the  cylinder  in  this  plane  is 
sufficient  to  generate  the  results  for  the  entire  cylinder  (Ref.  27). 

From  the  list  of  input  parameters,  the  fluid  sphere  is  idealized  into  FTAXI2 
MARTSAM  finite  elements  based  upon  the  number  of  elements  specified  for  the  sphere 

*  DATFFC:  from  data  generation  program  for  the  analysis  of  free-flooded  cylinders. 
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surface  and  the  cylinder  lengtli  and  thickness.  The  FTAXI2  element  is  a  compressible  fluid 
axisymmctric  clement  with  a  triangular  cross  section.  Each  element  has  six  nodal  points 
with  a  single  degree  of  freedom,  the  pressure,  associated  with  each  node.  Figure  6  shows  a 
typical  fluid-sphere  idealization.  The  triangles  are  calculated  by  means  of  a  triangularization 
routine  developed  for  a  contour  plotting  program  (Ref.  28).  Lists  of  the  (r,  z)  coordinates 
of  each  nodal  point  in  the  fluid  sphere  idealization  (COORDF)  and  the  fluid  material  para¬ 
meters  and  clement  assembly  data  (ELMNTF)  are  produced. 


The  finite-element  model  for  the  piezoelectric-ceramic  cylinder  is  generated  by 
modeling  the  top  half  of  one  stave  of  the  cylinder  with  PHEX20  MARTSAM  finite  elements. 
Again,  the  idealization  is  determined  from  the  list  of  input  parameters.  The  PHEX20 
clement  is  a  solid,  three-dimensional  element  with  20  nodal  points  for  each  clement.  Each 
nodal  point  has  three  displacement  degrees  of  freedom  (x,  y  and  z  in  a  rectangular  coordinate 
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Iljiurc  6.  C  toss-scclional  display  of  llic  finite-element  idealization  of  the  top  half  of 
the  acoustic  fluid  sphere  that  encloses  the  cylinder. 

system)  and  an  electrical-potential  degree  of  freedom  v.  Fig.  7  shows  a  typical  stave 
idealization.  Lists  of  the  (x,  y,  z)  coordinates  of  each  stave  nodal  point  (COORDS)  and  the 
piezoelectric-ceramic  material  parameters  and  stave  clement  assembly  data  (ELMNTS)  arc 
calculated. 

When  the  coordinate  data  COORDS  and  assembly  data  ELMNTS  arc  used  in  the 
M  ARTSAM  finite-clement  program,  the  mass  and  stiffness  matrices  arc  generated  for  the 
stave  model  in  a  rectangular  coordinate  system.  Bccuasc  an  axisymmctric  structural  repre¬ 
sentation  is  desired,  it  is  necessary  to  transform  the  equations  to  a  cylindrical  coordinate 
system.  The  DATFFC  program  generates  a  list  of  rotation  angles  for  each  nod.’i  point 
(BCLIST),  which  MARTSAM  can  then  use  to  perform  the  coordinate  system  transformation. 
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Hgurc  7.  Hnitc-clcr.ient  idealization  of  (he  top  half  of  one  stave  of  the  cylinder. 

Two  lists  arc  generated  that  will  provide  MARTSAM  with  the  necessary  information 
to  constrain  the  three-dimensional  stave  model  into  a  plane  (i.c.,  a  two-dimensional  axisym- 
mctric  model)  while  still  retaining  the  proper  electrical  properties  for  a  tangentially  polarized, 
staved  cylinder.  The  BAKADL  list  sets  equal  the  radial- and  vertical-displacement  degrees 
of  freedom  for  nodal  points  with  identical  r-  and  z-coordinatcs.  The  MPACKL  list  packs 
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out  the  appropriate  radial,  azimuthal,  and  vertical  degrees  of  freedom.  In  addition,  the  two 
lists  apply  the  appropriate  electrical-potential  boundary  conditions.  All  of  the  potential 
degrees  of  freedom  for  the  nodal  points  on  one  side  of  the  stave  are  set  equal  to  represent 
the  driven  foiled  surface,  and  all  of  the  potential  degrees  of  freedom  for  the  nodal  points  on 
the  opposite  side  of  the  stave  arc  set  equal  to  zero  to  represent  the  grounded  foiled  surface. 
Tlic  remaining  electrical-potential  degrees  of  freedom  are  left  free  and  arc  condensed  out 
later. 

The  program  DATFFC  also  generates  the  spherical  harmonic  transformation  matrix 
(LI  and  the  structure/acoustic  fluid  coupling  matrix  [Ti.  Finally,  lists  of  parameters  and 
matrix  sizes  are  produced  for  use  in  botii  the  second  and  third  programs. 

PROGRAM  MRTFFC* 

This  instruction  program  forms  and  assembles  the  mass  and  stiffness  matrices  of  the 
fluid  and  for  the  piezoelectric-ceramic  stave,  transforms  the  equations  for  the  stave  from 
rectangular  to  cylindrical  coordinates,  constrains  the  displacement  degrees  of  freedom  into 
an  axisymmctric  representation,  and  applies  the  appropriate  displacement  and  electrical- 
potential  boundary  conditions. 

Following  the  procedure  given  on  pages  10,  11,  the  electrical-potential  degrees  of 
freedom  arc  condensed  down  to  the  one  driving  potential.  Equations  (23)  and  (24)  are  then 
used  to  calculate  the  resonance  (short  circuit)  and  antiresonance  (open  circuit)  frequencies 
of  the  cylinder  in  air  (FREQST  and  FREQOT).  The  electrical  capacitance  (CAP)  in  air  is 
also  calculated  from  Eq.  (30). 

The  eigenvalues  and  eigenvectors  of  the  interior  problem  arc  calculated  from 
equation  (57),  and  the  KAUG  lowest  cigcnfrcquencics  (FREQY)  should  contain  at  least 
one  frequency  above  the  frequency  band  of  interest.  Augmenting  the  surface  equation  by 
too  many  degrees  of  freedom  docs  no  harm  except  to  slightly  increase  computational  costs, 
but  augmentation  by  too  few  equations  can  cause  singularity  problems. 

Finally,  the  matrices  that  comprise  Eqs.  (70)  and  (74)  are  assembled  for  use  in  the 
third  computer  program. 

PROGRAM  FRQFFC** 

This  program  requires  the  following  input  parameters  to  be  specified  by  the  user: 

IFFQ,  the  lowest  frequency  of  interest 
ILFQ,  the  highest  frequency'  of  interest 
IFQINC,  the  frequency  increment 

R(I),  distances  chosen  for  the  calculation  of  the  transmitting  voltage 
response  and/or  directivity  patterns 

*  MARTS  AM  program  for  the  analysis  of  free-flooded  cylinders. 

**  Frequency -sweep  computer  programs  for  the  analysis  of  free-flooded  cylinders. 


Using  Eq.  (74),  the  computer  program  FROFFC  calculates  the  real  and  imaginary 
parts  of  the  electrical  impedance  Z  for  each  frequency.  Also  at  each  frequency,  Eq.  (78)  is 
used  to  calculate  the  radiated  power,  and  then  the  solution  vector  is  normalized  to  Z  and 
the  transmitting  voltage  response  is  calculated  from  Eqs.  (79)  and  (80)  at  0  =  0  deg  and 
0  =  9  0  deg.  The  0-dcg  response  is  calculated  off  of  the  top  of  the  cylinder  along  its  symmetry 
axis,  while  the  90-deg  response  is  calculated  off  of  the  side  of  the  cylinder  on  its  symmetry 
plane  (see  Fig.  5). 

Vertical  directivity  or  beam  patterns  can  be  calculated  for  any  frequency  and  distance. 
Pressures  at  nodal  points  that  arc  interior  to  the  fluid-sphere  surface  and  the  structural  dis¬ 
placements  can  also  be  calculated  for  any  frequency.  From  this  data,  contour  plots  can  be 
made  of  the  magnitude  and  piiase  of  the  near-field  pressure  distribution. 


4  EXAMPLE  ANALYSES  AND  COMPARISONS  WITH  EXPERIMENTS 

To  assess  the  accuracy  of  predications  made  with  the  mathemat:  al  model  presented 
in  this  report,  three  existing  free-flooded  cylinders  were  analyzed.  These  cylinders  were 
chosen  because  they  have  quite  (iiffcrcnct  cliaractcristics,  and  experimental  measurements 
previously  made  of  the  electrical  impedance,  transmitting  voltage  response,  and  directivity 
patterns  were  readily  available. 

Cylinder  A  is  constructed  from  32  staves  of  barium  titanatc  (Ceramic  B  type) 
piezoelectric  ceramic.  To  waterproof  the  bare  ceramic  clement,  the  cylinder  was  fiberglass 
coated.  Tin's  was  done  to  keep  to  a  minimum  the  discrepancies  between  the  physical  trans¬ 
ducer  and  the  mathematical  model  of  the  transducer.  Except  for  the  waterproof  coating, 
no  other  additions  were  made  to  the  bare  ceramic  clement  before  testing.  The  agreement 
between  this  experiment  and  the  predicted  results  will  be  the  best  possible. 

Cylinder  B  is  a  60-stave  cylinder  constructed  from  a  type  of  lead  zirconatc-lcad 
titanatc  ceramic  (a  PZT-8  variety).  This  transducer  consists  of  the  ceramic  cylinder  wrapped 
with  fiberglass  to  provide  mechanical  bias  that  will  protect  the  ceramic  from  fracture  when 
driven  at  high  levels,  inner  and  outer  rubber  boots,  brass  end-rings,  and  castor  oil  surround¬ 
ing  the  cylinder  between  the  inner  and  outer  rubber  boots.  (Sec  Fig.  8  for  a  cross- 
sectional  sketch  of  a  cylinder  in  a  typicai  case.)  The  cylinder  can  now  be  easily  mounted, 
driven  at  high  levels,  and  is  electrically  insulated  from  the  water  by  the  rubber  boot  and 
castor  oil. 

Cylinder  C  is  a  30-stave  cylinder  constructed  from  another  type  of  lead  zirconate- 
lead  titanatc  ceramic  (a  PZT-4  variety).  This  transducer  again  was  tested  with  the  piezoelectric- 
ceramic  cylinder  wrapped  in  fiberglass  and  installed  in  a  case  similar  to  the  one  described  for 
Cylinder  B. 
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CYLINDER  A 


Cylinder  A  was  constructed  only  for  the  purpose  of  comparison  with  theoretical 
predictions.  The  32-stave  barium  titanate  (ceramic  B  type)  cylinder  is  not  wound  with 
fiberglass  filaments  nor  does  it  liavc  a  protective  case  or  any  provision  for  mounting.  Com¬ 
parisons  between  experiments  conducted  with  this  cylinder  and  theoretical  predictions 
sliould  yield  the  best  possible  agreement. 

Tabic  1  gives  some  of  the  parameters  used  as  inputs  to  the  computer  programs.  As 
was  pointed  out  in  Ref.  1 1 ,  material  parameters  that  arc  listed  in  handbooks  (Ref.  29)  are 
often  very  inaccurate.  Ceramic  varies  from  manufacturer  to  manufacturer,  and  differences 
in  material  parameters  as  great  as  20  percent  are  not  uncommon  in  supposedly  identical 
types  of  ceramic.  In  an  attempt  to  produce  material  parameters  which  more  accurately 
describe  the  ceramic  out  of  which  the  cylinder  under  consideration  was  constructed,  a  simple 
perturbation  technique  was  used  to  fit  the  breathing  resonance,  breathing  antiresonance, 
and  tire  electrical  capacitance  of  the  cylinder  in  air.  In  general,  this  technique  works  well. 
Tiie  book  and  perturbed  values  for  the  piezoelectric-ceramic  material  parameters  are  listed 
in  Table  2.  Table  3  gives  the  measured  values  for  the  electrical  capacitance  and  the  first 
three  resonance  and  antiresonance  frequencies  of  the  cylinder  in  air  and  also  the  predicted 
values  on  the  basis  of  both  the  book  parameter  values  and  the  perturbed  parameter  values. 


Table  i .  Input  parameters  for  the  three  example  free-flooded  cylinders. 


Parameter 

Cylinder  A 

Cylinder  B 

Cylinder  C 

RI  (m) 

0.1641 

0.09843 

0.08573 

RO(m) 

0.1927 

0.I15I 

0.09287 

RLNTII(m) 

0.1272 

0.08636 

0.07938 

NSTAVE 

32 

60 

30 

SMALLA  (m) 

0.22 

0.14 

0.12 

RHO  (kg/m3) 

1000 

1000 

1000 

RCV  (m/scc) 

1500 

1500 

1500 

NPSR 

9 

9 

9 

NLR 

2 

2 

2 

NTR 

1 

I 

1 

NSPIIH 

10 

10 

10 

KAUG 

10 

10 

10 

The  measured  and  predicted  real  and  imaginary  parts  of  the  in-water  electrical 
impedance  arc  shown  in  Rg.  9.  The  small  discrepancies  between  the  computed  and  the 
measured  responses  are  attributed  to  inaccurate  materia!  para.^-ters  and  limitations  in  the 
model  (inclusion  of  no  material  losses  and  the  inability  to  model  glue  joints,  for  example). 


perturbed  pie/.oeicctric-oeramic  material  parameter 


Tlie  measured  and  predicted  transmitting  voltage  responses  are  shown  in  Fig.  10. 
Hie  measurements  and  calculations  were  made  at  a  distance  of  1  m  off  of  the  side  (6  =  90- 
deg)  of  the  cylinder.  The  agreement  is  again  excellent. 
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Figure  10.  Transmitting  voltage  response  of  Cylinder  A  as  a  function  of  frequency. 

Solid  line  is  for  experiment  and  dashed  line  for  theory. 

Figures  1 1  and  1 2  display  the  measured  and  predicted  vertical  directivity  or  beam 
patterns  at  3,000  and  7,000  Hz,  respectively. 

The  maximum  source  level  and  maximum  radiated  pow  er  can  be  calculated  by 
assuming  the  maximum  permissible  electrical  field  will  occur  for  an  applied  voltage  of 
5V/0.001  in.  (2,000  V/cm,  a  conservative  estimate).  For  Cylinder  A,  the  minimum  distance 
between  adjacent  electrodes  is  3.2  cm  and  the  maximum  voltage  is  therefore  about  6,400  V. 
At  3,000  Hz,  the  maximum  source  level  is  2 1 3  dB  re  1  pPa  at  1  m,  while  at  7,000  Hz,  the 
maximum  source  level  is  210  dB  re  1  pPa  at  1  m.  The  minimum  source  level  for  the  cylinder 
driven  at  6,400  V  over  the  frequency  band  from  1 .5  to  10.0  kHz  occurs  at  about  5,000  Hz 
and  is  190  dB  re  1  n Pa  at  1  m.  The  maximum  radiated  power  at  3,000  Hz  is  6  kW,  and  the 
maximum  radiated  power  at  7,000  Hz  is  1 .7  kW. 

CYLINDER  B 

Cylinder  B  is  constructed  from  60  staves  of  lead  zirconate-lead  titanate  (PZT-8  type) 
tangentially  polarized  piezoelectric  ceramic.  The  in-air  measurements  of  the  lesonance  and 
antiresonance  frequencies  and  electrical  capacitance  were  performed  on  the  bare  ceramic 
element,  but  the  in-water  experiments  were  conducted  with  the  cylinder  wrapped  in  fiberglass 
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and  mounted  in  the  case  shown  in  Fig.  8.  Table  1  gives  the  dimensions  of  the  ceramic 
cylinder,  Table  2  lists  book  and  perturbed  piezoelectric-ceramic  material  parameters  for  the 
PZT-8  like  ceramic  material,  and  Table  3  shows  the  measured  and  predicted  in-air  response 
of  the  cylinder.  Again,  the  use  of  perturbed  material  parameters  produces  good  agreement 
between  theory  and  experiment. 


The  measured  and  predicted  real  and  imaginary  parts  of  the  electrical  impedance  are 
shown  in  Fig.  1 3.  The  agreement  between  experiment  and  theory  for  this  cylinder  is  not  as 
good  as  it  was  for  Cylinder  A.  The  presence  of  the  case  appears  to  have  two  major  effects 
on  the  response  of  the  transducer.  First,  the  case  ir  creases  the  length  of  the  cylinder,  giving 
the  cavity  a  greater  volume;  therefore,  the  first  cavity  resonance  occurs  at  a  lower  frequency 
than  predicted.  Second,  the  rubber  boots  and  castor  oil  introduce  larger  losses  into  the 
system  than  were  present  in  Cylinder  A.  Because  the  mathematical  formulation  used  for 
these  predictions  includes  no  structural  damping,  at  frequencies  where  losses  m  the  rubber 
and  castor  oil  are  appreciable,  the  predicted  impedance  does  not  compare  very  well  with 
the  measured  impedance. 


Figure  14  shows  the  measured  and  predicted  transmitting  voltage  responses.  Again, 
the  lack  of  damping  in  the  mathematical  model  and  the  increase  in  cavity  size  cause  some 
discrepancies.  In  addition,  it  appears  that  minor  resonances  of  the  case  or  mounting  hardware 
have  been  introduced  into  the  system. 


Figures  15  and  16  show  the  measured  and  predicted  vertical  directivity  patterns  at 
the  frequencies  of  4,500  and  6,000  Hz,  respectively.  The  agreement  of  4,500  Hz  is  excellent 
and  the  agreement  at  6,000  Hz  is  good.  The  differences  in  geometry  and  damping  between 
the  bare  ceramic  cylinder  and  the  transducer  as  tested  appear  to  produce  the  greatest  effect 
off  of  the  top  of  the  cylinder  (for  small  angles  of  0). 


The  maximum  source  level  and  maximum  radiated  power  can  be  calculated  in  the 
same  manner  used  for  Cylinder  A.  Assuming  a  maximum  voltage  of  2,000  V  the  maximum 
source  level  is  approximately  206  dB  re  1  pPa  at  1  m  at  4,000  Hz  and  193  dB  re  1  juPa  at 
1  m  at  8,000  Hz.  The  maximum  radiated  power  at  3,000  Hz  is  5  kW,  at  4,000  Hz  is  1  kW, 
and  at  8,000  Hz  is  about  Vi  kW. 


CYLINDER  C 


Cylinder  C  is  a  30-stave  cylinder  constructed  from  lead  zirconato.-lead  titanate 
(PZT-4  type)  piezoelectric  ceramic  with  the  dimensions  given  in  Table  1 .  The  in-air  measure¬ 
ments  of  the  resonance  and  antiresonance  frequencies  and  electrical  capacitance  were  per¬ 
formed  on  a  bare  ceramic  element,  and  quite  interesting  results  were  found.  From  Table  3, 
the  measured  in-air  response  of  the  transducer  shows  that  the  bending  resonance  and  anti¬ 
resonance  frequencies  lie  between  the  breathing  resonance  frequency  and  the  breathing 
antiresonance  frequency.  This  was  verified  by  using  an  optical  probe  to  measure  the 
displacement  distribution  along  the  length  of  the  cylinder  at  the  frequencies  of  interest. 

The  calculated  in-air  response  made  with  book  material  parameter  values  also  predicted  this 
phenomenon  but  not  at  the  correct  frequencies.  Performing  the  same  simple  three-parameter 
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Figure  15.  Vertical  directivity  pattern  for  Cylinder  Bat  4,500  Hz.  Solid  line  is  for  experiment 
and  dashed  line  for  theory  Pressures  are  plotted  in  dB  re  1  /iPa  for  1  V  input  at  1  m. 
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Figure  16,  Vertical  directivity  pattern  for  Cylinder  B  at  6,000  Hz.  Solid  line  is  for  experiment 
and  dashed  line  for  theory.  Pressures  arc  plotted  in  dB  re  >  jUPa  for  1  V  input  at  1  m. 
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perturbation  calculation  that  was  performed  for  Cylinders  A  and  B  resulted  in  only  a 
partial  success.  The  calculated  resonance  and  antiresonance  frequencies  more  closely  fitted 
the  measured  frequencies,  but  the  bending  resonance  and  antircsonance  frequencies  no 
longer  lay  between  the  breathing  resonance  frequency  and  the  breathing  antiresonance 
frequency.  It  appears  that  cither  the  simple  perturbation  calculation  performed  here  is  not 
sufficient  to  fit  all  of  the  material  parameters  or  the  finite-element  model  of  the  cylinder 
lias  too  few  elements  to  predict  the  proper  bending.  However,  this  discrepancy  in  frequencies 
is  only  about  5  percent  and  should  not  affect  the  in-water  response  of  the  transducer.  Table  2 
i;-.es  the  book  and  perturbed  piezoelectric-ceramic  material  parameter  values. 

Figure  1 7  shows  tiie  measured  and  predicted  real  and  imaginary  parts  of  the  electrical 
impedance.  The  experiment  was  performed  with  the  ce»amic  cylinder  again  housed  in  a  case 
similar  to  Fig.  8.  The  predicted  impedance  closely  follows  that  measured  experimentally 
except  at  the  two  frequencies  where  again  damping  ."attens  the  measured  response. 

Tiie  measured  and  predicted  transmitting  voltage  response  for  Cylinder  C  is  shown 
in  Fig.  18.  Except  for  the  overestimation  of  the  level  at  the  cavity  resonance  frequency,  the 
two  curves  show  excellent  agreement. 

Figures  19,  20  and  21  display  the  measured  and  predicted  vertical  directivity  patterns 
at  3,000,  6,000  and  10,000  Hz,  respectively.  Just  as  for  the  previous  cylinder  housed  in  a 
case,  the  measured  and  predicted  patterns  do  not  agree  as  well  at  0  ~  0  deg,  as  they  do  at 
0  =  90  deg.  The  added  structure  and  the  presence  of  structural  damping  give  the  transducer 
a  response  that  is  slightly  altered  from  that  of  the  bare  ceramic  elements. 

The  maximum  source  level  and  maximum  radiated  power  are  found  by  assuming  a 
maximum  voltage  of  3,500  V.  The  maximum  source  level  is  therefore,  for  example,  203  dB 
re  1  pPa  at  1  m  at  3,500  Hz  and  falls  to  only  197  dB  re  1  pPa  at  1  m  at  8,000  Hz.  The 
maximum  radiated  power  is  7.4  kW  at  3,000  Hz  but  down  to  0.4  kW  at  5,750  Hz  and  to 
0.2  kW  at  8,000  Hz.  These  are  again  conservative  figures  since  acceptable  drive  levels  are 
often  6  or  8  V  per  0.001  in. 


5  CONCLUSIONS 

The  powerful  mathematical  modeling  techniques  that  have  been  presented  in  this 
report  can  be  used  to  calculate  the  response  of  any  axisymmetric  transducer  (or  transducer 
that  can  be  approximated  by  an  axisymmetric  representation  over  the  frequency  band  of 
interest)  when  immersed  in  an  infinite  acoustic  fluid  medium.  In  particular,  this  report 
describes  the  computer  programs  that  have  been  written  to  use  this  mathematical  formula¬ 
tion  in  the  analysis  of  tangentially  polarized  piezoelectric-ceramic,  free-flooded  cylinder 
transducers.  The  programs  are  designed  to  take  a  few  basic  parameters  and  produce  a  set 
of  graphs  that  display  the  complex  electrical  impedance,  transmitting  voltage  response,  and 
radiated  power  as  functions  of  frequency.  At  any  selected  frequency,  the  vertical  directivity 
pattern  (the  horizontal  directivity  pattern  is  always  circular  since  an  axisymmetric  response 
is  assumed)  and  displacement  distribution  on  the  cylinder  can  be  calculated. 
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Figure  18.  Transmitting  voltage  responae  of  Cylinder  Cm  function  of  frequency. 
Solid  line  I*  for  experiment  and  dashed  line  for  theory. 


Figure  20.  Vertical  directivity  pattern  for  Cylinder  C  at  6.000  llz.  Solid  line  is  for  experiment 
and  dashed  line  for  theory.  Prc/aires  arc  plotted  '»  *.i  re  i  #ll*a  for  I  V  input  at  I  m 


mathematical  modeling  difficulties  but  would  increase  computational  costs.  At  present 
computer  rates,  the  minimum  total  cost  to  operate  the  computer  programs  described  in  this 
report  is  approximately  S12S.  This  high  cost  for  the  analysis  of  each  proposed  free-flooded 
cylinder  limits  the  design  capabilities  of  the  model.  However,  as  a  final  check  or  a  last 
iteration  on  a  set  of  dimensions  chosen  by  some  simpler  and  more  economical  means,  the 
mathematical  modeling  techniques  described  in  this  report  offer  predictions  that  are  nearly 
as  accurate  as  experimental  measurements  and  are  much  cheaper  and  faster. 
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